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Abstract 

We consider the totally asymmetric simple exclusion process 
(TASEP) in discrete time with sequential update. The joint distri- 
bution of the positions of selected particles is expressed as a Fredholm 
determinant with a kernel defining a signed determinantal point pro- 
cess. We focus on periodic initial conditions where particles occupy 
cffl, d > 2. In the proper large time scaling limit, the fluctuations 
of particle positions are described by the Airyi process. Interpreted 
as a growth model, this confirms universality of fluctuations with flat 
initial conditions for a discrete set of slopes. 

1 Introduction 

The totally asymmetric simple exclusion process (TASEP) in discrete time 
consists of particles on Z with at most one particle at each site (exclusion 
principle). At each time step, particles jump to the neighboring right site 
with probability p £ (0,1), provided the target site is empty. There are 
mainly two types of update rules one can consider. One is called parallel 
update. It consists in first checking for all particles if they can jump (i.e., if 
their right neighbor is empty) and then, simultaneously and independently, 
these particles jump to the right each with probability p. The second update 
rule, the one we actually analyze in this paper, is called sequential update. 
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In this case, at each time step particles are processed sequentially from right 
to left, i.e., starting from right to left, if the site on the right of a particle 
is empty, then this particle jumps there with probability p. This update 
rule allows to shift blocks of particles to the right during one time-step. 
Equivalently, it can be regarded as a parallel update rule for holes, which 
under the constraint of keeping their order, jump to the left k steps with 
probability proportional to p k . 

On a macroscopic level the particle density, u(x,t), evolves deterministi- 
cally according to the Burgers equation dtu + d x (u(l — u)/(p~ l — u)) = 0. 
It is therefore natural to focus on fluctuations, which turn out to have un- 
expected features. For example, the fluctuations of particle positions live on 
a t 1 / 3 scale and the limiting distribution depends on the type of the initial 
configuration. 

In this paper we consider the positions of particles as observables. If 
one starts with step initial conditions, i.e., particles occupying Z_, then 
the macroscopic particle density has a region of linear decay [14]. There 
the fluctuations are governed by the GUE Tracy- Widom distribution [17], 
see [7]. Moreover, in the limit of large time t, the joint distributions of particle 
positions are governed by the Airy process [12] (called Airy2 process in the 
following). The second situation which has been analyzed is the stationary 
initial condition [6]. These results have been obtained for the continuous time 
and/or parallel update discrete time TASEP, but can be easily extended to 
the sequencial update case. 

A third class of initial conditions consists in deterministic and periodic 
initial configurations. In [15] and [1] the continuous time TASEP has been 
studied for the particular case of alternating initial configuration, i.e., initially 
with particles on TL. In this case the one-particle fluctuations are described 
by the GOE Tracy- Widom distribution. Joint distributions are given by 
a Fredholm determinant and the kernel converges pointwise to the one of 
the Airyi process. In the present paper we extend these results in several 
directions. Firstly, we look at the fully discretized version of the TASEP 
with sequencial update. Secondly, our result holds for a wider set of periodic 
initial configurations and, thirdly, the control on the kernel is stronger and we 
prove convergence of joint distribution functions. The corresponding results 
for continuous time can be obtained along the same lines. 

The Airyi process, Ai, is a marginal of the signed determinantal point 
process with extended kernel K Fl , in the same way as the Airy 2 process is 
related to the extended Airy kernel. Explicitly, we set B (x,y) = Ai(x + y) 
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and let A denote the one-dimensional Laplacian. Then 



K Fl (u 1 ,s l ;u 2 ,s 2 ) = -(e^- u ^ A )( Sl ,s 2 )l(u 2 > Ul ) + ( e - WlA J B e^ A )(s 1 , s 2 ). 

(1.1) 

More explicitly, as shown in Appendix A of [1], one has 

K Fl (ui,si;w 2 ,s 2 ) = 7= :exp - — l(u 2 > «i) 

+Ai(si + s 2 + (w 2 - Mi) 2 ) exp {{u 2 - uO(si + s 2 ) + ^{u 2 - u x fj . (1.2) 

In particular, the one-point distribution of the Airyi process is given in terms 
of the GOE Tracy- Widom distribution [18], as shown in [5], 

P(A(0) <s) = F x (2s). (1.3) 

The process A\ was first described by Sasamoto [15]. The starting point 
is a determinantal formula for the probability distribution of the continu- 
ous time TASEP with finitely many particles, first discovered by Schiitz [16]. 
Sasamoto found a clever reformulation making possible the large time asymp- 
totic analysis for the positions of particles. The details of the derivation are 
presented in [1]. There we show how the signed determinantal point process 
arises by applying a result on the L-ensembles defined in [2] . 

For the discrete-time sequential update TASEP, a determinantal formula 
of the same type as in [16] was derived in [13] via the Bethe Ansatz. Since this 
is the starting point of our analysis we present here an elementary derivation. 
Very recently, a determinantal formula for parallel update was obtained [10]. 
Whether or not in this case a similar approach can be used for asymptotic 
analysis has still to be investigated. 

By universality, one expects that the same limit process appears for any 
deterministic and periodic initial configurations, not only for particles ini- 
tially on 2Z. In this paper we consider initial conditions with one particle 
every d sites. The fluctuations of the position of a particle live on the scale 
t 1 / 3 , while the positions of two particles are non-trivially correlated (on the 
t 1//3 scale) if they are of order t 2 ^ 3 apart. In Theorem 12.51 we prove that, 
properly rescaled, the joint distributions of particle positions converge to the 
joint distributions of the Airyi process. That is, let Xk{t) be the position at 
time t of the particles starting at — dk, then for some constants //, a, £i,£ 2 
depending only on p and d, one has 

- t:: — -rjr, > Ai[u), as t — > oo, (1.4) 

at i 6 
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Figure 1: The growth model associated to the TASEP. If a particle jumps to 
the right, then the surface growth vertically as indicated by the dashed line. 



in the sense of finite dimensional distributions. 

Reformulation of the result. The TASEP can be interpreted as a stochas- 
tic growth model of an interface, which turns out to belong to the Kardar- 
Parisi-Zhang (KPZ) universality class introduced in [8]. Given a realization 
of the particle process, one defines for each time t a continuous height func- 
tion, x I— > h t (x), as follows. The height at is given by h t (0) = 2N t , where N t 
is the number of particles which jumped from site to site 1 during the time 
interval [0, i). If at time t a particle is at site x, then h t decreases linearly 
with slope —1 in the interval [x,x+ 1], and increases with slope +1 if the site 
is empty. Thus, if x is an integer, one has always \h t (x + 1) — h t (x)\ = 1. A 
particle jumping at time t from x to x + 1 corresponds to increasing h t (x + l) 
by 2. The surface obtained in this way grows vertically under the TASEP 
dynamics as indicated in Figure [TJ 

The step-initial conditions, corresponding to h (x) = \x\, generate a 
curved macroscopic shape and the limit process is the Airy2 process, first de- 
scribed in another growth model, the polynuclear growth model [12]. The al- 
ternating initial conditions, particles starting from 2Z, correspond to growth 
on a flat and horizontal substrate, h Q (x) = (1 + (— l) x )/2, x G Z. In this case 
the macroscopic limit shape is a constant function and the fluctuations are 
governed by the Airy! process. By universality one expects that the result 
are unchanged if the initial substrate is flat but tilted. Theorem 12.51 confirms 
universality for a discrete set of non-zero slopes (corresponds to d > 3). 

The TASEP can also be interpreted as a directed percolation model, 
where "flat" corresponds to the point-to-line, while "curved" to the point- 
to-point setting, see e.g. [4,7,11]. Finally, it has been conjectured [1,3] that 
the evolution of the largest eigenvalue of GOE Dyson's Brownian Motion of 
random matrices is governed by the Airyi process too. This conjecture is 
however not based on KPZ universality. 
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New features. In the analysis some new interesting features appear. The 
kernel can be expressed in a double integral representation and, as stated in 
Theorem 12.11 one of the contour integrals circumscribes d — 1 simple poles, 
which are the roots of a polynomial. There is only one simple case, d = 2, 
which is equivalent to flat and horizontal initial conditions. In this case, 
there is only one pole and the asymptotic analysis reduces to the one of a 
single integral on the complex plane. An explicit solution of the polynomial 
equation does not exist for arbitrary d > 3, thus we had to employ several 
new ideas to circumvent the problem. Similar situations might appear in 
other problems like point-to-line directed percolation and all sorts of growth 
models on a flat substrate. 

Outline. The remainder of the paper is organized as follows. In Section [2] 
we define the model precisely and state the main result of this work. In 
Section [3] we give a new derivation of the determinantal formula of [13] and 
then apply our previous work [1] to the discrete-time TASEP. In Section H] 
we perform the necessary orthogonalization and obtain the finite time kernel. 
The asymptotic analysis is the content of Section [51 In Appendix [A] we deal 
with the trace-class problem of the kernel Kp 1 . In Appendix [B] we describe 
an alternative, more constructive, way of obtaining the orthogonalizazion. 

Acknowledgment 

A. Borodin was partially supported by the NSF grant DMS-0402047 and the 
CRDF grant RIM1-2622-ST-04. P.L. Ferrari thanks M. Loss for sketching 
the proof that the one-point kernel Bq is trace-class during his visit to Munich 
in 2005. 

2 Model and results 

The model analyzed in this paper is the discrete-time TASEP with sequential 
update on Z. At any given time t, every site j G Z can be occupied at most 
by one particle. Thus a configuration of the TASEP can be described by 
V — {Vjyj Z|7/j G {0,1}}. r]j is called the occupation variable of site j, 
which is defined by rjj = 1 if site j is occupied and rjj = if site j is empty. 

Let r](t) be a TASEP configuration at time t. Then the configuration at 
time t + 1 is obtained by the following dynamics. Starting from right to left, 
a particle jumps to the neighboring site with probability p G (0, 1) provided 
this site is empty. Since the update is sequential from right to left, during a 
time step a block of consecutive particles can jump. For example, if at time 



5 



• — • 1 — - • 1 • - — I • — • ► 

1 — p p(l — p) p 2 

Figure 2: Example of the transition weights with sequential update. Above 
is the configuration at time t. Below are the possible outcomes at time t + 1 
together with the transition weights. 



t we have r) (t) = r]i(t) = 1 and rj 2 (t) = 0, then the particle at can jump to 
1 provided the particle at 1 jumps to 2, see Figure[21 

Let us start at time t — with particles occupying the sublattice dIL. The 
observable we focus at is the joint distribution of particle positions, which 
can be expressed as a Fredholm determinant as stated in Theorem 12.11 

Theorem 2.1. Let <j(l) < cr(2) < . . . < a(m) be the indices of the m 
particles starting at time t — from positions 2^(0) = — d(a(k) — 1). Denote 
by Xk{t) their positions at time t. Their joint distribution is given by 

m 

P ( fl {^(*)(*) - = det ( 1 ~ XaKXa)e2({*(l),...M m )}xZ), (2- 1 ) 

k=l 

where Xa{o~{k),x) = t(x < a^) and the kernel K is defined as follows. 

For fixed v, letuo(v), . . . ,Ud-i{v) be the roots of the equation u(l+u) d ~ l = 
v(l + v) d ~ l , with uq{v) = v the trivial solution. Then the kernel for the d- 
spaced initial configuration is 

K(n 1 ,x 1 ;n 2 ,x 2 ) = -\ Xl %2 ]) + K (n 1 ,x 1 ;n 2 ,x 2 ) (2.2) 

V^ 2 - ni - I J 

with 

Ko(ni,Xi;n 2 ,x 2 ) (2.3) 

1 I dv y^ l + dv (l+pUi^Yi-Uiiv^il + v)™' 2 
2iri J ro V j^l + dui{v) (1 + jw)*(-v) n2 (l + Ui(v)) Xl+ni+1 ' 

where T is any simple loop anticlockwise oriented around the pole at v = 
(without other poles being inside T ). 

For alternating initial configurations, d — 2, the kernel is particularly 
simple. 
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Corollary 2.2. The kernel for the 2-periodic configuration is given by 



K(n 1 ,xi;n 2 ,x 2 ) 



+ 




xi+ni+ri2— 1 



X2+ni+Tl2—2 



1 — p — pv 



1 + pv 



) 



(2.4) 



Proof. Apply (12.31) with U\(v) 



1-v. 



□ 



In Theorem 12.51 we prove the convergence of the fluctuations of particle 
positions to the Airy! process. Thus we need to define this process and the 
scaling limit. 

Definition 2.3 (The Airy! process). Let B (x,y) = Ai(x + y), A the one- 
dimensional Laplacian, and the kernel Kp 1 defined by 



K Fl (u 1:Sl ;u 2 ,s 2 ) = -(e^~ u ^ A )(s u s 2 )l(u 2 > Ul ) + (e-^ A B e^ A )(s u s 2 ). 



The Airyi process A\ is the process with m-point joint distributions at 
u\ < u 2 < . . . < u m given by the Fredholm determinant 



where Xs(u k ,x) = l(x > s k ). 

Remark 2.4. The kernel Kp 1 is not trace-class on L 2 ({ui, . . . , u m } x R), 
because the diffusion part appearing for u 2 > u\ makes it not even Hilbert- 
Schmidt. However, as shown in Appendix [XJ there exists a conjugate 
operator, one with a kernel leading to the same Fredholm expansion of 
det(l — XsKf-lXs), that is trace-class on L 2 ({ui, . . . , u m } x R). Thus the 
Fredholm determinant in i\2.6\\ regarded as its Fredholm expansion series is 
well defined. 

We focus on the region around the origin. The fluctuations of the particle 
positions are of order t 1 ^ 3 and non-trivial correlations occur for particles at a 
distance of order t 2//3 , as expected from the KPZ scaling exponents. The mean 
density of particles is 1/d and the probability of jumping to the neighboring 
site is p, provided the position is available. With sequential update, a particle 
can jump to the right even if the target site is occupied, provided that the 
blocking particle itself jumps during the same time step. Therefore in the 



(2.5) 



rn 



P( f]{Mu k ) < = det(l - XsK ¥l x s )m{ 



Ui,...,u m }xH); 



(2.6) 



fc=i 
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stationary state the probability that the right position is available is given 
by 

The average speed of particles is p(d — l)/(d — p), which means that at time 
t particles with index close to tp(d— l)/((d — p)d) will be close to the origin. 
This motivates the following scaling. 

Scaling limit. Define the constants 

(2(l- P )p)V3 W _i)) 2 /3 2 

^ ' ''--"d^r (2 - 8) 



We consider particles with index n(u, t) given by 

n(u, t) 



d(d — p) d 



(2.9) 



which typically at time t is at a position close to jjut 2 ^ 3 . Recall that the 
initial condition is £fc(0) = — d(k — 1), k £ Z. Then the rescaled process is 
given by 



u i-> X t (u) = _^ 1/3 (2.10) 

and converges to the Airy! process as follows. 

Theorem 2.5. Let X t be the rescaled process as in li2.1U\) . Then, in the limit 
of large time t, it converges to the Airyi process Ai, 



\imX t = A 1 , (2.11) 



t— >oo 



in the sense of finite- dimensional distributions (given by 112.6)) ). 



3 Signed determinant al point process 

We first consider a system with a finite number N of particles. The <i-periodic 
configuration is then obtained by the proper iV — > oo limit. We start at 
time t — with N particles at positions yjv < • • • < 2/2 < Vi and study the 
probability that at time t these particles are at positions x^ < ■ ■ ■ < %i < 
Denote this transition probability by 

G t {xi,...,x N ) = P{(x N , . . . ^^^Kun, . . . ,yi,0)). (3.1) 
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Lemma 3.1 (See [13]). The transition probability has a determinantal form 



G t({x}) = det(F i ^ j (x N+ i- i - yN+i-j,t))i<i :j <N 



with 



F n {x,t) = {l-p) 



(-1)' 
2m 



dw 

I 0,1 W 



p 



1 — p 



-W 



(3.2) 



(3.3) 



where T^i is any simple loop around and 1 oriented anticlockwise. 

Notational remark: with the notation J Fg f(z)dz with S a set of points 
we mean that the integral is taken over any simple loop, oriented anticlock- 
wise, enclosing the set S, and such that all the poles of / inside T$ belong 
to S. Equivalently, J r f(z)dz = 27ri ^] 2g5 Res(/, z). 

This result is already contained in [13], where the derivation used the 
Bethe-ansatz method. Here we present a different derivation. 

Proof of Lemma \3.1\ To obtain (13.21) we first write the master equation for 
the iV particles, 



G t+ i(xjv,...,xi) = E 

6i,...,6jv6{0,1} 



N 



P 



1 — p 



rb N 



G t (x N -b N ,...,x 1 -b 1 ) 
(3.4) 

with the boundary conditions due to the exclusion constraint 

Gt(%N, ■ ■ ■ , Xk+l, Xk+l + 1, %k-l, ■ ■ ■ , %l) =G t (X]y, . . . , Xfc+i, Xfc+i, Xfr-l, . . . , X\), 

(3.5) 

for k = 1, . . . , N - 1, cp. with [16]. 

We assume that a formula of the form (13. 2\\ holds and prove that there 
exists a family of functions {F n } satisfying the Lemma. As an abbreviation 
we set F^(x, t) = F n (xN+i-j— yN+i-j, t) during the proof. Inserting (13.21) into 
the master equation (13.41) and using the mult i- linearity of the determinant, 
we get 



Gt+li^N, ■ ■ ■ , %l) 
6i,...,6 Jv e{0,l} 



p 



p 



b!+---+b N 

det [F^xx+x-i - b N+1 _i,t 



l<i,j<N 



det 



E 

bN+l-i=0 



P 



P 



bjV+l-i 



(3.6) 



l<i,j<N 



det [(1 -pjF/^zjv+i-ijO +pi 7 T_ j (a;iv+i_i - 1,£ 
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l<i,j<N 



This is equal to det [i^ -(zjv+i-i, t + 1 



l<i,3<N 



if the functions F n 's satisfy 



F n (x,t + 1) 



[l-p)F n (x,t)+pF n (x-l,t). 



(3.7) 



We also have to take into account the boundary conditions. The two sides 
of f)3.5p are N x N determinants with all the rows identical except for the 
(N + 1 — fc)th one. Combining the two determinants, (13.51) reads 



= det 



F 3 

r N+l-k 



F 3 N _ k _j(x k+1 ,t) 
-j(xk+i, t) - Fj f+1 _ k _ j (x k+1 + 1, t) 



(3i 



l<j<N 

The two lines explicitly written are linearly dependent whenever the functions 
F n satisfy 

F n _ x (x, t) = c(F n (x, t) - F n (x + I, £)) (3.9) 

for some c. Here we choose c = 1. Any other choice of c corresponds to 
the replacement F n by c~ n F n , which, however, keep the determinant in (13. 2p 
unchanged. 

The functions F n are determined by the two relations (13. 7ft and (13.91) . as 
well as the initial condition 



Xi, 



8„ 



VN,%N 



(3.10) 



F (x,t) is already determined by one-particle configurations. In fact, in 
this case, G t (x) = P(x(t) = x\x(0) = y) = Fq(x — y,t). Therefore 



F {x-y,t) = (1-pY 



P 



1 — p 



x-y 



t 

x-y 



(3-11) 



This result is consistent with (13.71) and (13.101) . Denote by A the discrete 
derivative Af(x) = f(x + 1) — f(x). Then H 3 . 9 1) implies 



(3.12) 



F_ n (x,t) = (-l) n {A n F )(x,t). 
F has the following integral representation, 

F (x,t) = (l-py±<[ —U + T^-wY^, (3.13) 

27T1 J To W \ 1 — p J W x 

where T is any simple loop around oriented anticlockwise. Therefore to 
obtain F_ n we simply apply 



A'' 



w 



w) 

n+x 



(3.14) 
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t(-l) n I dw ^ p V (1- w) 



Thus, for n > 0, 

F_ n (z, t) = (1 - pf^^- & — (1 + ^— w) x ^ ~ J . (3.15) 
V ' ; V yJ 2m J Fo w \ l-p J w x+n V ; 

In this case, there is no pole at w — 1, thus replacing T by r 01 leaves the 
result unchanged. 

For n > 0, F n is determined by the recurrence relation 

F n+1 (x,t) = J2 F n(y,t) (3-16) 

y>x 

together with the property that Fq(x, t) = for x large enough. To have the 
sum in (13.161) well defined we need \w\ > 1, i.e., the integration path includes 
1, which is a pole of the integrand for n > 1. Thus in order for (13. 16ft to be 
satisfied for all n, we need to take into account the poles both at and 1. □ 

The analogue of Theorem 12.11 for a finite number of particles is the fol- 
lowing Proposition, which is a consequence of Lemma 3.4 of [1]. 

Proposition 3.2. Let a(l) < a(2) < ... < a{m) be the indices of m out 
of the N particles initially at y a (i), • • • , ya(m) ■ The joint distributions of their 
positions x a ^) (t) is given by 

m 

P ( fl { X °(k){t) < a k } ) = det(l - XaKXa)e2({a(l),...,a(m)}xZ) (3.17) 
k=l 

where Xa{o~{k),x) = t(x > a^). K is the extended kernel with entries 

"2-1 

K(n 1 ,x 1 ;n 2 ,x 2 ) = -<j>^\x 1 ,x 2 ) + ^ Kl^+iM^M ( 3 -18) 

i=0 

where 

( j>^\x 1 ,x 2 )=( Xl - X2 - 1 \ (3.19) 
\n 2 -n x -\) 

The functions n > 1, % 6 Z ; are defined by 

«T(») = (1- „)< ^(l + ^Vii^ (3.20) 

where the path Tq in the definition of is any simple loop, anticlockwise 
oriented, which includes the pole at w = but not the one at w = 1. The 
functions $™, i — 0, . . . , n— 1, n > 1, are polynomials of degree at most n — 1 
uniquely defined by 

Y,*i{x)*?(x) = 6 i)S , j = 0,...,n-l. (3.21) 

xew, 
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Proof. The proof is an application of Lemma 3.4 of [1], with 



<t> n (x?,xp 1 ) = l{x?>x] +1 ), n = l,...,iV-l, (3.22) 

and 

= (-l^-'F-N+iix -y h t), i = l,...,N. (3.23) 

The above functions Fi are defined by an integral enclosing w = only, since 
w = 1 is not a pole for i < 0. With the definition (I3.20p above we have the 
composition rule 

tf*V£i- j )(x) = Vl_ j (x). (3.24) 

In our setting, if we sum up all the variables {x™, 1 < m < n,l < j < m}, we 
obtain a Vandermonde determinant in the variables x™. Thus the space V n of 
Lemma 3.4 in [1] is generated by {1, x, ... , x n ~ 1 } and the are polynomials 
of order at most n — 1. A simple computation using ( 13.201) leads to 

= | ' " ' (3.25) 

which, together with (13.211) leads to $q(x) = 1 = 0n-i(oo, a;). Then the 
Proposition 13.21 follows from Lemma 3.4 of [1]. □ 



4 Orthogonalization 

Consider the case where the particles are initially regularly spaced as follows 

Vi = -d(i-l), i = l,...,N. (4.1) 

for d > 2. For any fixed time t, we first obtain the kernel for a fixed number 
of particles, N. Then we take the N — > oo limit and focus on a point far 
enough to the left where the particles do not feel the fact that there is a 
rightmost particle. 

Lemma 4.1. The functions *&k(x) and $^(x) have the following integral 
representations . Let z = x + d(n — 1), then 

and 

*»(x) = tHl dv 1 + dv f4 3) 

fcl J 27Ti jfp V (1+JW)* CKl+^) d - 1 ) fe 

where Tq is any anticlockwise simple loop enclosing only the pole at 0. 
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Proof. We have 

= {-l) k F- k {x-y n - k ,t) = (-l) k F. k (z-dk,t). (4.4) 

Thus (13.31) leads directly to (14.21) . Next we prove that (14.31) satisfies the 
orthogonality relation (13.2ip . Since ^ k (x) = for x < —d(n — 1), i.e. z < 0, 
we have 

j2n(^m(^)) = <f dw(i+ P (w - i)Y{(w - i)^- 1 )' 



z>0 



X ■ 



l) j f dv (1 + dv) ^ (v + 1) 



z-1 



2vri 7 ro u (1 + pvY(v(v + l)^ 1 )^ ^ io 

provided that the integration domain satisfies |1 + v\ < \w\. The last sum 
gives 

y <" + 1 >'" = I . (4.6) 

w z+1 (w - (1 +v))(l + v) v ' 

(14.51) has a simple pole at w — 1 + v. Therefore the integration over w leads 
to 



(—] ) k +j C 1 _|_ fin) 

= V— f dv ^r^vM l + v ) ) j - ( 4 - 7 ) 



^_!)fc+i /• i + dv 

V k \x{z) )<!>■■ [x{z)) = - 

2>0 

The final step is a change of variable. Let u = v(l + v ) d ~ l . Then 

du = (l + v) d ~ 2 (l + dv)dv (4.8) 

and the integral is again around 0. Thus 

.l)k+j r | 

z>0 



E ViMzWM*)) = <f L_ = 5 , ifc . (4.9) 



□ 

There is also a more constructive way of doing the orthogonalization using 
Krawtchouk orthogonal polynomials. In fact, as shown in Appendix [HI the 
vI/W's can be written as linear combinations of Krawtchouk polynomials on 
the interval [0, t + d(N — 1)] and parameter p. What one has to do is to invert 
a certain N x N matrix with entries depending on the initial conditions. We 
were able to do it in the case d = 2 and obtain the formula of Lemma 14. 11 
Once the form of the $'s is obtained for d — 2, it is easy to find the proper 
ansatz in the case d > 2. 
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Proof of Theorem \2.1[ We apply Proposition 13.21 with the orthogonalization 
of Lemma 14.11 With the change of variable u = w — 1 in (14 .2p the function 
has an integral representation with a pole at u = — 1, namely 

= {1 ^ )z+1 (l+pu) t (u(l + u) d - 1 ) k (4.10) 

with z = x + d(n — 1). The main part of the kernel, denoted by 
K (ni,x 1 ;n 2 ,x 2 ), reads 

n Jl} r_1^i-«a /• a + c/f)(l+w) Z2+d ~ 2 



k=0 ^ > JV ° 

J , (1 + pu)\u(l + u)^ 1 )™ 1 1 



(1 +pv) t (v(l + v) d - 1 ) n 2 



'l + u) Zl+1 u(l + uV- 1 - v(l + V 



with 2j = Xi + d(rii — 1) provided that the paths of integration for u and v 
satisfy (a) \u{l + w) d_1 | > \v(l + i;)^ -1 1 and that (b) u — is not inside 
the contour r_j. The expression (14.1 ip is obtained as follows from (13.181) . 
First we exchange the finite sum and the integrals. Then we extend the 
sum over k to — oo. The expression remains unchanged if the two above 
conditions are satisfied. In fact, for k < — 1, the pole at v = vanishes 
because of the extra term v ~ k and by (a) the series is absolutely summable. 
Condition (b) reflects the fact that in the definition of (I3.20p the integration 
path includes only the pole at 0. A possible choice of the integration paths 
is r_i = {u, \ l+u\ = 1/2} and To a contour with \v\ small enough to satisfy 
(a). 

The kernel for the ci-spaced initial configuration is obtained from the finite 
system by choosing the x« such that the kernel becomes independent of the 
fact that we have only a finite number of particles. In other words, we choose 
Xi such that the kernel is invariant with respect to the shift of the positions 
by — d and the particle numbers by 1. This is achieved when z% < (d — l)rii 
because u = — 1 is not anymore a pole. Then only the d — 1 simple poles 
inside T_x contribute. 

When we evaluate the integral over u we first integrate out the d — 1 
non-trivial poles, that is, the zeros of 

R(u, v) = w(l + uf- 1 - v(l + v) d -\ (4.12) 

Let Ui(v), . . . , Ud-i(v) be the solutions of R(u,v) = different from the 
trivial one Ud{v) = v. Then the kernel for the <i-spaced initial configuration 
is as given in (I2.2p - (I2.3I) . A simple case is d — 2, where the only non-trivial 
solution of R(u, v ) = is u = — 1 — v . □ 
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-L L £ t 2 / 3 oo 
Figure 3: The different regimes of (s\,s 2 ) for which the kernel is analyzed. 



5 Asymptotic analysis 

Proof of Theorem 12.51 



To prove Theorem 12.51 we consider a conjugate kernel to K, i.e., a kernel 
which gives the same correlation functions. It is given by 



K co «i{n 1 ,x 1 ;n 2 ,x 2 ) 



K(m, x\\ n 2 , x 2 ) 



(5.1; 



d- 1 



X2+dri2 — (x\+dn{) 



(d-1) 



d-1 



We also use the notation K^ 011 - 1 for the conjugate of the second part of the 
kernel K, denoted by K in (12. 3p . Notice that the kernel is invariant under 
simultaneous shifts Xi — > X{ + dS and rii — > Hi + S , S £ Z. For the scaling 
limit we have to consider, compare with (12.81) and ( \2.9L 



11; 



d-l 



^ - 



d(d — p) d(d — 1) 



t + K 2 



£/ 3 



(5.2) 



From now on we fix r 1; . . . , r m £ R. The dependence on these constants 
is not indicated in the following propositions, since it is irrelevant for the 
proof of Theorem 12.51 

To prove convergence of the Fredholm determinant we need to control 
the behavior of the kernel (15 .ip as a function of Si and s 2 . We do it in three 
steps corresponding to the three regions indicated in Figure El In Propo- 
sition 15.11 we prove uniform convergence on a bounded set, by controlling 
the finite-t deviations from the asymptotic term. We thus have the con- 
trol for (si,s 2 ) £ [— L,L] 2 for any fixed L. Then, in Proposition 15.31 we 
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obtain a bound for (si, s 2 ) G [— L, et 2//3 ] 2 \ [— L, L] 2 , for L and t large enough 
and e small enough. Finally, in Proposition 15.41 we obtain a bound for 
(si, s 2 ) G [— L, oo) 2 \ [— L, et 2 / 3 ] 2 , for any fixed e > and £ large. 

Proposition 5.1 (Uniform convergence on compact sets). For fixed r\ and 
r 2 , the extended kernel has the following limit. Let us fix any L > 0. Then, 
with Xi,ni defined as in $5.2\) . the kernel converges uniformly for (si,s 2 ) G 
[-L, L} 2 as 

lim K con i(n u x l ;n 2 ,x 2 )Kt 1/3 

t—*oo 

exp f- [ S2 ~ Sl)2 ) l(r 2 > n) (5.3) 



/ 2 
+Ai(si + s 2 + (r 2 - n) ) exp (Jr 2 - ri)(si + s 2 ) + -(r 2 - ri 



Proposition EH] leads to a uniform bound on the kernel as in Corollary l5.21 
This is obtained using the super-exponential decay of the Airy function, 

2 3 /2 

namely, for any a > 0, |Ai(x)| < e~3 x < Const a e _aa; for some Const a 
depending on a. 

Corollary 5.2. For any /ixed L > tnere exists a to — ^o(-^) > and a 
Const ^ independent of to s.t. fort > t the bound 

11^(711, xi;Ti2,a:2)Kt 1/3 1 < Const L (5.4) 

holds for all si, s 2 G [— L, L] . 

Proposition 5.3 (Moderate deviations). For any large enough L there exist 
an eo = Sq{L) > and a to — t (L) > such that for any positive s < and 
t > to the estimate 

\K con \n l ,x 1 ;n 2 ,x 2 )Kt 1/3 \ < e - (si+Sa) (5.5) 
holds for (si,s 2 ) G [-L, et 2 / 3 ] 2 \ [-L, L} 2 . 

Proposition 5.4 (Large deviations). Let c = (d — l)p/(d — p), e > as in 
Proposition \5.3\ (small enough), and s~i = KSit~ 2 ^ 3 . 

(1) If §i G [e,c] or s 2 G [e,c] we have the bound 

\K^{n 1 ,x 1] n 2 ,x 2 )\<e-^ + ^ (5.6) 

for t large enough. 

(2) If §i G (c, oo) or s 2 G (c, oo), then the kernel is zero for t large enough. 
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The last ingredient is an estimate on the binomial part of the kernel. 
Proposition 5.5. For any s\, s 2 G R and r 2 — r\ > fixed, the bound 



Kt 



1/3 



,d-ly V(tZ-l)*- 1 / \ n 2-ni-l, 

< Const ie - |s2 - Sl1 (5.7) 

/ioWs /or t large enough and Consti independent of t. 
With these results we can now prove Theorem 12.51 

Proof of Theorem \2.5[ From Corollary 15 .21 Propositions 15.31 and !5.41 it follows 
that there exist a to > and a Const2 independent of to such that 

\K c on \n 1 ,x 1 -n 2 ,x 2 )Kt 1 / 3 \ < Const 2 e- (siV0+S2V0) (5.8) 

for any t > to. Here a V b = max{a, b}. 

The joint distributions of the rescaled process X t defined in ( 12.101) are 
given by a Fredholm determinant with series 



P[f]{X t (u k )<s k }^ 

k=l 

/_i\n m r " 

!" Yl d yi ...dy n ]J l(y k < ^u lk t 2 ' 3 - KS ih t x ' 



k=l 

i \n m r n 

1/3^ 

n>0 ii,...,i n =l 17 fc=l 

x det (iT on >(u ifc ,t), b/*];n(« ill t) l N)) 1<fcKn - (5.9) 



By the change of variables o k = (y k — fxui k t 2 ^ 3 ) / (— nt 1 / 3 ) and a conjugation 
we obtain 

m „ n 

(JE9])= J2 £ /do- 1 ...d ( 7 n J]l( ( 7 fe > Si J (5.10) 



n>0 ii,...,i n =l fc=l 



x^det [^Mu^t), N;^K,t), + 



Kfc.Kn 



The term t^^t^ is the new conjugation, which does not change the deter- 
minant. 

On i^Q° nj we use the bound (15.81) . Whenever i k > ii, we have an additional 
term coming from the binomial part of the kernel. This contribution is 
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bounded as in Proposition 15.51 Therefore the (k, I) coefficient in the last 
determinant in (I5.10p is bounded, for t large enough, by 



Const 3 (e-l^l + e-^ v0 +^ v °))£g|^, i k > i h 
Const 3 e-^ v0 +^ v °)^J, i k <i h 



(5.1i; 



with Const 3 independent of t. We also use the bounds, for i, j e {1, . . . , m}, 
{ Pr^je-^ < - i 2 £± ^ eH^I < -^Const 4 , j > i, (5.12) 

(l + irp 1 + y (i + y ) 1 + 2/ 

and 

(l + ^ 2 ) t -xvo-,vo < Const e -(*-h/)/2. (5.13) 
These bounds applied to (15. lip leads to 

(IM < Consta { £ + ^ (5.14) 

for some Constg independent of £. 

(I5.14p implies that in the determinant in (15.101) we can single out a prod- 
uct f3^ =1 max{(l + al)~ 1 ,e~ CTk ^ 2 }. Then the entries of the determinant are 
bounded by Const6, so that the whole integrand is bounded by 

1 n 

-ConstgV*/ 2 JJmax{(l + ^)- 1 ,e- CTfc/2 }l(a fe > —L) (5.15) 
n ' fc=i 

where the factor n n l 2 is the Hadamard bound (the absolute value of a de- 
terminant of a n x n matrix with entries of absolute value not exceeding 1 
is bounded by n n ^ 2 ). Therefore, for t large enough and some Const7 > 
independent of t, 

E n/2 
— -m n Const" < 00 (5.16) 
n\ 

n>0 

because n! grows like (n/e) ra , which is much stronger than n n l 2 and any term 
exponential in n. Therefore for t large enough the integrand can be bounded 
by a t-independent integrable function. By applying dominated convergence, 
we can exchange the t — > 00 limit with the sums and integrals. The pointwise 
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convergence comes from Proposition 15.11 thus 

m 

lim v(C\{X t {u k ) < s k } 

k=l 

_^ (_]\n m ^ r " 

= ~ S / d0_1 • • - dcPn II det (- K Fi( U <fc> <7 fc; u *H t7 »))l<fe,i<n 

n>0 ' h,...,i n =l ak>Si k k=l 
m 

= p(f]Mi(« fc ) <«*})■ ( 5 - 17 ) 

fc=i 

Thus X 4 converges to as i — > oo in the sense of finite-dimensional distri- 
butions. □ 

Proof of the uniform convergence on bounded sets 

Proof of Proposition \5. 11 Let us consider the scaling (15.21) . First we determine 
the contribution coming from the binomial term of ( 12.2ft . i.e., 

/ j \ xa+dna-fci+dni) / , d x ni-na / -,\ 

-(z4t) ((^r (::::;)• 



Set 

2k 2 



d(d- 1) 



(r 2 -ri)t 2/3 -l,6 = /<s 2 -si)t 1/3 + (d-l), and £ = b/a. (5.19) 



The binomial term is ( ^„ ) and x! = V 27rx exp(x ln(x) — x)(l + (9(x x )). 
Therefore we have 

( 15~T8D = 1 / +0(rl/3) e -fc 2 /(2^-i)q) (1 + ^(e)). ( 5.20) 
V47r(r 2 - ri) 

By reinserting the expressions for a and b we obtain 

ASHED = . =i = ex P M?~ ai) n (i + 0(e) + 0(r 1/3 )) (5.21) 

with an error uniform for Si,s 2 G [—L,L\. For any fixed r 2 > r 1; the error 
O(e) is also of order 0(t~ 1 ^). Thus taking the limit t — > oo we get the first 
term of the limit kernel Kp 1 , see (12. 5p . 

Next, we consider the main part of the kernel, i.e., (12. 3ft . It is not difficult 
to see (by looking at the image of any closed simple loop around u — — 1 
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under the map w(l+M) a! ~ 1 ) that the image of {ui(v), . . . , Ud-\{v)} form closed 
loops around — 1 if v is a closed loop around 0. The leading term of the kernel 
is 

1 d_1 f 

— V* dvexp{til)(v)-ti/j(ui(v))) (5.22) 
2711 i=i -^0 

= -—E—\ n ( v ) _ln(l+pu). (5.23) 
a — p 



with 



We can rewrite 

ReMjfv)) = -rr^—]n(\v\\l/p + v\^ p - 1 ) -lnp. (5.24) 
a — p 

In Proposition 15.91 we prove that there exists a path T of finite length such 
that Re(ip(v)) — Ke(ip(ui(v))) < whenever Ui(v) ^ v. Moreover, Proposi- 
tion [523 says that Re(V>(v)) —Re(ijj(ui(v))) < except at the point v = —1/d 
for the solution U\(y) (the one for which u±(— 1/d) = = i>). There- 

fore the contributions of « 2 , • • • , are bounded by O(e~ Uit )F(0) for some 
Vi > and 

((^wj d^i- (5 - 25) 



The contribution of Mi can be estimated by integrating only on \v + l/d| < S 





for some 5 > 0. The error will be exponentially small in t. Thus, let Tq 



denote To restricted to a 5- neighborhood of —1/d. Then, 

K {ni, xi; n 2 , x 2 )/tt 1/3 = O( e -^)F(0) (5.26) 
+ KtV3 / 1 + ^ ( l +P u ^ v ))\- u ^ v )) ni ( l + v ) X2+n2 ~ 2 



2id J r s l + dux(v) (l+pv) t (-v) n *(l + u 1 (v)) x i+ n i+ 1 

for some fi = /i(5) > 0. 

From now on we denote u — u\ since the other solutions of R(u, v) = 
do not appear anymore. Close to the critical point we can compute the value 
of u by the series 



1 , 1. 2d(d-2), 1,2 4d 2 (d-2)% i 3 i 



^) = ___( B+ _ ) __^ (B+ _ ) ^__i_ f( „ + _)» +0((l , + _)- ) . 

(5.27) 

Let z = v + 1/d and let the image of Tq be t/j, where \z\ < 5. Thus we can 
use Taylor expansion and the last term of ( 15. 26ft becomes 

^- J d^exp (tf (z) + t^h{z) + t 2 ^f 2 (z) + h{z)) (5.28) 
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with 

d 3 K 3 
3(7^1) 3 



/ 2 (z) = -K(s 2 -s 1 )H(d-l)/d)-K-^—(s 1 + s 2 )z + 0(z 2 ), 

a — 1 

/ 3 (z) = ln(d/(d-l)) + 0(z). (5.29) 

The errors are uniform for Sj G [— L, L] (the s- dependence in the error terms 
is only in the j 2 term). Let fi(z) be the expression fi(z) without the error 
terms. Define F(z) = exp(tf (z) + t 2/3 /i(^) + t 2 ^f 2 (z) + h{z)) and F(z) 
similarly. Then 

Kt 1 / 3 r Kt 1 / 3 r ~ Kt 1 / 3 r 

^ I d2F ' 2 > = ^ I + I <**W - F W ). (5.30) 

To estimate the last term we use the inequality | — 1 1 < el x l|a;|. Thus 

Kt 1 / 3 



< 



-2vri 

Kt 1 ' 3 



dz(F(z) - F(z)) 



ns 



(5.31) 



[ dz\F( Z ) \ e °(**t+zW*+zW' + z) ( zH + z 3 t 2/3 + ^1/3 + ^ 



= Ktl/3 f ^| e t/oW(l+Xl)+t 2/3 /lW(l+X2)+t 1/3 /2W(l+X3)| 

x 0(zH + /t 2 / 3 + zh 1 ' 3 + z) 

for some xi, X2, X3 which can be made as small as desired by choosing 5 small 
enough. We take as integration path rjs = {e~ 17rsgn< - w ^ 3 \w\,w G [—6, 5]} which 
is close to the steepest descent path when w — > 0. 

At the integration boundaries, w = ±5, the leading term is 

thus the integral remains bounded as t — > 00 and, by the change of variable 
wt 1 / 3 = q we see that the error term is (9(£ _1 / 3 )F(0). 

The final step is to compute 3^- J dzF(z). Extending 5 to 00 we collect 
only an error of order O(e~ tlt )F(0) with < fi ~ 5 3 . This leads to the 
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integration on the path r]^ = {e mssa ( w ^ 3 \w\,w G R}. Thus 

nt 1 ' 3 f ~, , _ N 2/rf 1 / 3 e- i ^ 3 f , / d 3 K 3 



dsF(*) = F(0) — / dwexp ( - — — — u > d t 



- 27ri Aoc " ~ 27ri V 3(d-l) 

x exp ( "'[^y e-^V^) exp (-/^(si + ^e-^tirt 1 ^ . 

The change of variable g = wne~ m ^ 3 d/(d — 1) gives 

^! / cb%) = G(0)— *— [ dge <?3/3 e (r2 - ri)g2 e- (si+S2)<? (5.32) 

with G(0) = F(0)(d— l)/d. By considering the conjugate kernel K™ n> instead 
of K , the term G(0) cancels. 

Finally we use an Airy function representation 

— *— f dve v3/3+av2+bv = Ai(a 2 - b) exp(2a 3 /3 - ab) (5.33) 
— 27ri L 

•J '/oo 

to obtain the final result 

lim ^lK (x uni ; x 2 , n 2 ) = Ai( Sl + s 2 + (r 2 - Ti yy^-r 1 )>/3+(s 1 +. a )( n -r 1 ) 
t^oo G(U) 

(5.34) 

uniformly for si, s 2 G [— L, L]. □ 

The goal of the following sequence of lemmas is Proposition 15.91 used in 
the proof of Proposition 15. II above. 

Lemma 5.6. Define the path 

7! = {-1 + sin(0(d - l))e*/ sin(0d), G [0, 7r/d)}. (5.35) 

On 7 1; t>(l + i;) a!_1 G R. Lei U{, i — 1, . . . ,d, be the d solutions of the equation 

R(u, v) = u(l + uf' 1 -v(l + vf- 1 = 0. (5.36) 

Then, for all v G 71 and Wj G" {f , we /iaf e 

(^(l/p + Mi)^- 1 ) ^ ^(l/p + v)^ 1 ! . (5.37) 

Lemma 5.7. On 71 \ {—1/d} the solutions u{v) of R(u,v) = are simple 
zeros. 
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Lemma 5.8. Assume that a, 6 G K satisfy either 
(1) b < a < or 

f^J a < 6 and a(a — 1) < 6(6 — 1). Then, for /j, > 1, we have the strict 
inequality 

/i a -/ < 0- l)(a-6). (5.38) 

Proposition 5.9. Let u\(v) be the solution of R(u,v) = such that 
u\(—l/d) = —l/d, and U2(v), . . . , Ud-i(v ) the other non-trivial d — 2 so- 
lutions. Then, there exists a path T encircling the origin, passing through 
v = —l/d such that 

Re(^(f)) - Re(V>(ui(«))) < (5.39) 
except for the solution Ui(v) at the point v = —l/d. 

Proof of Lemma 15.61 For d = 2 nothing has to be shown because there are 
only two solutions, u = v and u — v. Thus we consider d > 3. Since for 
v G 71, v(l + v) d ~ l G R, we have to prove that the system of equations 

Hli + wl^- 1 = Hii + wi 6 ^ 1 , . 

\u\\l/p + u\ d l' p ~ l = \v\\l/p + v\ d / p - 1 ^ ' 

has only the trivial solutions u = v and u = v. If (15. 40 ft has a solution u, then 
by symmetry u is also a solution. Thus in the rest of the proof we restrict 
ourselves to the upper-half plane. 

(I5.40p is an equation involving only distances of u and v from the points 
0,-1, and —1/p. Let us choose any v and set 

a — \v\, b = \1 + v\, c=\l/p + v\. (5-41) 

Similarly we choose ati/« and set 

d=\u\, b=\l + u\, c=\l/p + u\. (5.42) 

Then (E30D writes 

~ab d - 1 = ab d -\ ~a~c d/p ~ l = ac d/p -\ (5.43) 

that is, 



aji, 



6/ Ai 1/(d - 1) , c = c /^ /{d - p) . (5.44) 



For n — 1, we have u = v. Thus consider /i ^ 1. Given a, 6 such that the 
circle centered at of radius a intersects the circle centered at —1 of radius 6, 
the position of v is uniquely determined, see Figure HI Thus this determine 
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-1/p -1 
Figure 4: Geometric representation of ( 15.40ft . 



c too, i.e., there exists a function / to be determined such that c 2 = f(a, b). 
Define the function 

Fr>, a, b) ee /f>, ft/zx 1 /^!)) _ /( a , b)/iFlM. (5.45) 

Then since also c 2 = f(a,b), the scaling relations (15.441) imply that the 
Lemma will be proven once we show that F(fi, a,b) = only at fi = 1 for 
a, 6 such that t> G 71. 

To determine the function / we use elementary trigonometry. Let 7r — j3 
be the argument of v. Then 

b 2 = a 2 + 1 - 2acos(/3), 

c 2 = a 2 + p~ 2 - 2ap~ 1 cos(/3) . (5.46) 



From this it follows 



c =f(a,b) = — 5 1 a. (5.47) 

p A p p 



v G 71, is parametrized by an angle and we have 



sin(0) ft = sin(0(rf - 1)) ^ 



sin (0c?) ' sin(0d) 
Therefore 

,2 



FQi, a,b) = ^ (1 - /i"*/^)) + 1 ~ PSin J^ - M 2 ) 
p 2 v 7 p sin (0a) 

+ sin 2 (0(rf - l)) _ 2/(d _ 1} _ (5 4Q) 

psin (0a) 
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Define the function 



= ^(i-/.-*^>) + ^0«-*«^>-^ 

" ,_1)2 (^-W-U - (5.50) 



po? 2 

Using the properties 

1) s "a /f.l and sm /tll"^ 1 ^ are positive and increasing function in 6 for 
<f>e [0,7r/d), 

2) ^- 2 P/( d -p) - /i 2 > for /i G (0, 1) and - /i 2 < for p > 1, 

3) p- 2 /^- 1 ) - ^p/id-P) > for G (0, 1) and p- 2 /^" 1 ) - ^vl^-v) < 
for /i > 1, 

we obtain that 

F(fi,a,b) >G(fi), for ^6(0,1), 

F(fi,a,b) < G(fi), for /x > 1. (5.51) 

Therefore we have to prove that G(/i) > for p G (0, 1) and G(fi) < for 
H > 1. This follows from the fact that 

1) G(l) = (trivial verification) and 

2) G"(aO < for all \i ^ 1 (to be proven below), 
is given by 



Consider first [l > 1. We can rewrite 
#(//) = (d - 1) (^P/^-P) - p- 2 /^ 1 )) + (1 - p) ^-ZP/V-P) - ^) . (5.53) 
By Lemma [5.81 if follows, for \i > 1, that 

AM = aT 2 ^ - < - ^ (a* - 1). (5.54) 

Moreover, for d > 3, by Lemma [5.81 we have that 

-2p/(*-„) _ ^2 < __J^_ (;U _ i) (5.55) 
a — p 
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holds for /x > 1. Using these bounds we obtain H{p) < for /x > 1, from 
which G'(/x) < for /x > 1 too. 

Now consider the case /x G (0, 1). We get 

fiH'(fi)/2 = {^ d -^ - +(l-p) (^/^ - /x 2 ) > (5.56) 

because both terms are strictly positive for /x G (0, 1). Since H(l) = 0, we 
then get H(p) < for /x G (0, 1). Thus Gr'(/x) < for /x G (0, 1) too, and this 
finishes the proof of Lemma 15.61 □ 

Proof of Lemma \5. 7[ On 71 we have t> (1 + t> G R. By symmetry if u(v) is 
a solution of R(u, v) = 0, then also u(v) is a solution of the same equation. 
Thus we have to check that u = v is a simple zero for v G 71 \ {— ■ We 
have -R(w, v) = Ylk=i( v ~ u k( v )) ■ Hu = v would be a double solution, then the 
Taylor expansion of R(u, v) at u — v would have the form J2k>i a k{v){v — u) k 
with ax(v) = 0. However, explicit computations gives 

R(u, v) = -v(l + v) d - 1 (- + — ^ J (v - u) + a 2 (v)(v - uf + ■ ■ ■ . (5.57) 

\V V + 1 J 

We see that the only v G 71 such that ai(f) = is v = — 1/d. Thus 
Lemma 15.71 is proven. □ 

Proof of Lemma \5.8[ Consider the function /(/x) = /x a — /x b . Then its second 
derivative is 

f"(n) = (a(a - l)/x a - 6(6 - l)/x 6 )/x~ 2 . (5.58) 

For 6 > a and /x > 1 we then have f"(fi) < (a(a — 1) — b(b — l))/i a-2 < 
as soon as a(a — 1) — b{b — 1) < 0. In this case /(/x) is strictly concave for 
/x > 1 and case (2) of Lemma [5.81 is proven. 

For b < a < 0, /(/x) is not anymore concave for all /x > 1. Let us compute 
the zero of f'(fJ-), /xi, and the zero of /"(jtx), /X2. We get 

„, = exp > 1, (5.59) 

* = ex P (— J ex P ( J355^ J > 

This means that /(/x) is strictly concave as long as it increases and at /x = 1 
it is still increasing. Thus the bound for case (1) follows. □ 

Proof of Proposition \5.iA We first define r close to infinity and the critical 
point. Secondly, using Lemma [5.61 we prove that we can complete the path 
satisfying Proposition 15.91 The way we do it is illustrated in Figure [51 
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Figure 5: Illustration of the path r (the bold line) used in the asymptotic 
analysis. The dashed line is the 71. On the right the structure close to the 
critical point is shown. 



First we consider the infinity, i.e., v = pe lip for p ^> 1, and (p G 
(—7i/d,ir/d). Then the d — 1 non-trivial solutions Uk{v), k = 1, . . . ,d — 1, 
(different from the u = v one) of R(u, v) = can be expanded in power of 
p~ x . Let us set 

u k = pe^e- [2 ^ d (l + a k /p + 0(p' 2 )). (5.60) 
Simple calculations lead to 

a k = d lZlh _ (_ 1) 2(^)(d-i)/dN e -V ) (5 61) 

d ' 

from which 

Re(ip(v))-Re(ip(u k (v))) = - ^~ P ^ Re(e^-e^e- i2 ^ k/d )+0(p- 2 ). (5.62) 

pp(d — p) 

Since for k = 1, . . . , d - 1 and ip G {-ir/d, ir/d), Re(e iv - e^e" 12 ^'^) > 0, 
it follows that for p large enough, Re(ip(v)) — R< e (V'(' li fc(' t '))) < 0- This strict 
inequality holds also for (p = ±7r/ d for k = 2, . . . , d — 1. Therefore we can 
choose a part of To to be parametrized by v = pe llf for all <p G (—n/d, n/d). 

Next we look close to the critical point. The above strict inequality holds 
also for (p = ±ir/d for k — 2, . . . , d — 1, thus also on the line 71. By continuity 
and Lemma Re(^(f )) — Re(ip(uk(v))) does not change its sign along 71, 
thus we can follow the path 71 to come back to the critical point, with the 
result 

Re(V>(-l/d)) -Re{ip{u k {-l/d))) < 0, k — 2,...,d—l. (5.63) 
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By continuity, Re(?/>(t>)) — Ke(ip(uk(v ))) < in a small enough neighborhood 
of —1/d. Thus any path leaving from — 1/d is good, at least locally, for the 
d— 2 solutions. Next, we consider the solution u±(v) close to the critical point. 
It is the only non-trivial solution that u\{— 1/d) = —1/d (see Lemma 15.71) . 
By Taylor expansion we obtain (15.271) . which yields 

ij>(v) - HMv)) = 3( f^ 3 1)3 ^ 3 + 0(z 4 ), z = v + 1/d, (5.64) 

see (15.291) . Thus in a neighborhood of the critical point there are three lines 
along which Re(i/j(v)) — ~Re(ip(ui(v))) = 0. These lines are locally given by 

v = - + i x + o(x 2 ), v = -- + e in/e x + 0{x 2 ), v = -- + e i5n/6 x + 0{x 2 ), 
d d d 

(5.65) 

for x real, see Figure Therefore we can choose the path T close to the 
critical point to be any path of the form 

v = -- + e w x + 0(x 2 ), 9 e (-7r/2,-7r/6) U (n/6,n/2). (5.66) 

The steepest descent path leaves the critical point with an angle ±7r/3. 

The final step is to see that we can join the part of the path To close 
to the critical point and the one far away (at large enough distance p from 
the origin) by going close to 71. Assume that we can not do that. Then 
somewhere along the way we must hit a point where Re(ip(v)) = Re(ip(ui(v )) 
with Ui(v ) 7^ v . Let us call points v with this property "bad" . Thus, we know 
that there are bad points at arbitrarily small distances to 71, and also these 
points have bounded absolute values because the neighborhood of infinity is 
completely controlled. 

Therefore, the sequence of bad points necessarily has a limit point on 71 
which is not —1/d (because near —1/d everything is controlled, too). But 
because of Lemma 15.61 this implies that bad points near 71 can come only 
from the root u\(v) which is equal to v on 71. Thus we have to prove that 
the root U\ does not cause bad points. For this, the neighborhood of 71 away 
from the critical point can be parametrized as 

.W) = (-i + »)r« (5.67) 

V sm(a<^) / 

For ip G (0, tt/gT), 1 ^> 9 > means that the point has increased real and 
decreases imaginary part with respect to 9 = 0. We can compute explicitly 
U\{v) in series of 9, with the result 

Ul (<p, 9) = v(-<p, 9)(1 + a9 + 0(9 2 )) (5.68) 
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with 

1 + l+v(-^,0)d V ; 

Let v(ip, 0) = x + h/. Then explicit computations lead to 

PS/ 



^(ReV(^)-ReVK^))) 



6»=o (d — p)((l + x) 2 + y 2 )((l + pa;) 2 + p 2 y 2 ) 
x ((p + 2pd - 3d){y 2 + x 2 ) + (3p - 2 - d) + (4d + 2pd - 2 - 4d)x) 

The denominator is always positive. Consider y > (i.e., away from the 
critical point). Then p + 2pd — 3d — p(l — p) + 3(p — l)d < for p e (0, 1) 
and d > 2. Thus the term in ?/ 2 is always strictly negative. To analyze the 
contribution in x, we set x = — 1/d + z, and on 7x \ { — l/d} we have z > 0. 
As a function of z, the remainder of the numerator becomes 

(d - l) 2 (d - p) 2(d - l)(p(d - 1) - 2d(l - p)) 
I 2 d Z 

-{p{d- 1) + 3d(l -p))z 2 < 0. (5.70) 

Therefore, in a right-neighborhood of 71, Re(ip(v )) — Re(^(wi(u ))) < 0. Thus 
tti(i>) does not generate bad points in a right-neighborhood of 71 (excluding 
71). □ 

In the above construction of the contour T we ignored the requirement 
that the contour is not allowed to contain any additional poles, see Theo- 
rem EH1 These additional poles appear as Ui(v) = — l/d. In particular, if 
u(l + u)^ 1 is real, then also v (1 + v ) d ~ 1 is real. This happens on the d — 1 
branches, one being 7! from (15.351) . the others are d — 2 branches originate 
at —1 and going to infinity in the directions e ±lk7T ^ d , k = 2, . . . , d — 1. These 
branches do not intersect with 71 because of Lemma I5.7[ thus do not inter- 
sect the interior of Tq. Therefore the path To fulfills all the requirements of 
Theorem 12.11 

Proof of the bound for moderate deviations 

Proof of Proposition [3T31 In this proof we set Sj = KSit~ 2 / 3 . The e > is still 
to be chosen. We can set it as small as desired, but of course independent 
of t. For e > small enough, the path To can be chosen to be equal to the 
path used in the proof of Proposition 15.91 and illustrated in Figure El except 
for a deformation close to the critical point. 
We set 

p / 1 + pv \ 

i>e,i{ v ) = --, hy(v/Ui(v))-hi — +s 2 ln(l+w)-siln(l+Mi(f)). 

d-p \l+pUi{v)J 

(5.71) 
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First consider the e = situation. Let Ui(v) be the solutions of R(u, v) = 
as in Proposition 15.91 It follows from Proposition 15.91 that there exists a 
5 > such that Re(ip eil (v)) < -5 for all v G r \ {\v + l/d\ < Const 8 <5 2/3 } 
for some Constg > small enough. Next consider the situation we want to 
analyze, i.e., e > small. Since To remains on a bounded region, we can 
choose e with < e <C 5^ 3 such that Re(^ el (f)) < 5/2 for v as above. 
Similarly, for i = 2, . . . , d — 1, Re(i[) Ej i(v)) < 5/2 for all v G r , since the only 
solution u%(v) such that Re (ipo,i ( — 1/^)) = is Ui(v). 

Thus we choose Tq as in Proposition 15 . 91 for > c5 2 ^ 3 }, and modify 

only the path in the region {\v + l/d\ < c5 2 ^ 3 }. The contribution coming 
from the unmodified path is then bounded by 0(e~ 5t / 2 ) < 0(e~ 2e3/it ) < 
0(e~ 2e3/2t ). Now, note that sf 2 < e 3 / 2 t, thus the previous contribution is 
bounded by 0{e~^ +s ^ 12 ). 

Next we have to consider the neighborhood of v = —1/d and see how 
the positions of the critical points depends on e. We will then be able to 
choose the appropriate path locally. We have already computed u\(v) in the 
neighborhood of v — — 1/d in (15.271) . The doubly critical point at v — — 1/d 
is separates now into two critical points v± on the real axis. They are the 
stationary points of Re(ip £t i(v )) — Re(V ; £ ,i(— 1/d)) and can be computed in 
series of \/si + §2- Explicitly we have 

Re(^>))-ReOM-l/rf)) = ^J^^ ^- ^J 2 ^ z{\ + 0{z)) 

(5.72) 

with z = v + 1/d. From ( 15.721) it follows that 

v± = -i ± v^ + l^^l + 0{y/h+T 2 )). (5.73) 

Re(ip £t i(v)) is minimal at v + . Therefore we can choose the path To to pass 
by v + . Since we need just a bound, we can disregard the error term and set 
the path T to pass by 

«i = -i±V^Tl^# (574) 

The steepest descent path leaves u9 along the imaginary direction, see Fig- 
ure E], thus the final choice of To in the (^-neighborhood of —1/d is simply 

= {^+ + iy, \y\ < y/h + h ^a/^ } j° int witn Vi and Vi- 
We change the variable by setting v = —1/d+z, and denote z + = v+ + l/d. 

Then the contribution coming from the integral over To in the 5-neighborhood 

of —1/d is 

dzexp(F t (z)) (5.75) 

V»7oVr;i 
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Figure 6: Local choice of the path r (the bold line) used in the asymptotic 
analysis for moderate deviations. 



where F t (z) = tf (z) + t 2 / 3 h{z) + t l ^f 2 {z) + f 3 (z), with 

/oW ^ ^Al^-^z^Oiz)), 

2, x 2 ((d-l) d - l \ , K 2 (r 2 -n)d 2 2 , 3 , 



f 2 (z) = -K(s 2 - Sl )ln((d-l)/d), 

f 3 (z) = \n(d/(d-l)) + 0(z). (5.76) 

This is like (I5.29P except that now the are in the leading term f since 
they can be very large, namely of order 0(t 2 ^). 

On 770, z = y/Si + s 2 (1 + i^) with w G (— v^? \/3)- It is not difficult 
to check that on r/ 



Re(*/o(*)) < -^s 1 + s 2 f 2 (l + 0( v fhTT 2 )), 



Re(t 2/3 /iW) < ^ 2/8 A(Q) + ( si + s 2 ) {l + Q(y/h + T 2 )), 

MM*)) < fsW + Oiy/h + h). (5.77) 

For L large enough, the linear term in S1+S2 is controlled by the — |(si+S2) 3 ^ 2 
term. Thus for e small enough, the contribution (15.751) coming from r/o can 
be bounded as 

\ Vo \Kt l ^e F ^e-^ +s ^ /2 . (5.78) 

But I rjo I = 2t~ 1 ^^/s\ + S2 ^~F~^ ; therefore the t 1//3 factor simplifies and the 
contribution ( 15. 75ft coming from r/o is bounded by Constge' Pt( ' * ) e _ 's( Sl+ ' S2 ) 3/2 for 
some finite Constg > 0. 
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The final step is to bound the contribution coming from the integration 
on rji (the same bound holds for the fji contribution). It is parametrized 
by z = e 17r ^ 3 y, y G [z+/2,5]. The term t 2 ^ 3 fi(z) is dominated, for large 
L, by the cubic term in z coming from the term tfo(z). In fact, while y 
increases the cubic term becomes more and more dominating with respect to 
the quadratic term in z and from the analysis on the domination occurs 
already at y — Z+/2. 

The linear term in z of fz{z) is dominated by the linear term in z coming 
from tfo(z), again for L large enough. Therefore for L ^> 1 and 5 < 1, 

where the last inequality follows from y > z + /2. Thus 

Kt 1 ' 3 



-2vri 



dze Ftiz 



m 



< Const 10 e Ft (°)e-i( Sl+S2 ) 3/2 / dyt 1 / 3 e- t ^^ y 



(5.80) 

The last integral equals a constant independent of t, and the t-dependent 
terms in the factor e^*- -* vanishes if we consider the conjugate kernel i^Q° nj 
of K , see (15. lft . In fact, 



e F t (0) 



d d ) \ d J d-l 

— l) d_1 \ ni_ ™ 2 / d _ ]_\ X2+dn 2 -(x 1 +dn 1 ) 



d d V d I d-l 



(5.81) 



Let the e and the t chosen above be denoted by e and to respectively. 
Then the decay is exponentially small in (sj + S2) 2 ^ 3 times a constant inde- 
pendent of L,e and t for < e < e and t > t . Thus, the whole result 
can be simply bounded by e~^ 1+S2 ^ for L large enough. Therefore the bound 
(1531) holds. □ 



Proof of large deviation bound 



Proof of Proposition 5.4 ■ In this proof the notation ~ means that the two 



expressions are equal up to a factor which is not exponentially large in t. 
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More precisely, we say that / ~ g if lim^oo | ln(/ / g) = 0. We have 

n 2 -l /^/_l\ :El+ni / rid \ Ul 

Let us denote a = a(k) = k/t. Then we have 

ri2 — 1 

Ko 0nj (ni,xi;n 2 ,x 2 ) ~ A Si)I 2 {a{k), s 2 ) (5.83) 



fc=0 



where 



ii(a,s) ~ — / duexp (tf a , s (u)), 

I 2 (a,s) ~ — / dwexp ( - tf a>8 (v)), (5.84) 
2vri /r 

and 

fa,s(u) = ln(l + pit) + a ln(— iz) — (d — 1) ln(l + u) 

+s(ln(l + u)-ln(l-l/d)). (5.85) 

We just have to find bounds on and J 2 such that their product is ex- 
ponentially small in t. Then, since the sum ( 15.831) includes 0(t) products, 
the result is obtained by determining the a G [0,p/(d — p)] which minimizes 
h(a, s)I 2 (a, s). 

The stationary points of f a>s (u), are by the Cauchy-Riemann equations 
also the critical points of Ke(f a>s (u)). Denote (5 — (d — l)(p/(d — p) — a), 
then we have to consider only a > 0, resp. s < /3, because the limit cases 
a = 0, resp. s = (3, correspond to a vanishing pole at u — 0, resp. u — — 1. 
Thus when a = 0, I 2 (a, s) ~ 0, and when s > (3, h(a, s) ~ 0. This is what 
happens in case (2), because if Sj > c = (d — l)p/ (d — p), then (3 < s for all 
a > 0. 

First we consider s = 0. Then 

df a o(u) pad , 

J ' 0K ' = — — + — = 5.86 

du 1 + pu u 1 + u 

has two solutions in (—1,0), 

1 a(d — p) 



of p(cv(d — p) + (1 — p)) 



. 1 a(d — p) . 

" ^ { ~ d' - P (a(d- P ) + (1- P )) .^ = (5 ' 87) 
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with equality for a = (1 — p)p/(d — p) 2 G (0,p/(d — p)) for any p G (0, 1) 
and d>2. In Lemma [5. 101 we prove that for s G (0,(3), u^ tS G (—1, — 1/cZ) is 
strictly decreasing in s, while u +yS G (— l/d, 0) and is strictly increasing in s. 
u +tS is the left-most maximum of f a ,s(u) for w G (—1,0). Thus, for s > Sq, 

/a,s(^+,s) > /a,s(w+,s ) > fa,s {u+, s ) (5.88) 

because / a , s (w+, so ) = / a , So (u+, So ) + (s - s )(ln(l + w+, So ) - ln(l - 
Similarly we o btain, for s < s Q , / QjSo (w_ So ) < f a , s (u- tS ). 

In Lemma 15.111 we prove that Ji(a,s) ~ e^"' 3 ^ 1 -' 3 ' (resp. 12(0;, s) ~ 
e -t/ a , s (w +>s )^ f rom w hich 

Ko 0nj (ni,xi;n 2 ,a;2) ~ exp (t max (/a,5i(«-,«i) - */o,«(«+,ia)))- ( 5 - 89 ) 

0<a<p/(c(— p) 

Then using (I5.88P we get 

Ko 0n \n 1 ,x 1 ;n 2 ,x 2 ) < exp ( t max (/ Q)£ (u_ )£ ) - /a, e (i; + ^)) ) 

\ 0<a<p/(d— p) / 

< ^ V [- 2 -{2e/^fh + tO{e 2 )). (5.90) 

For e — ■> the maximum is obtained at a = (1 — p)p/{d — p) 2 . For e > 
small it turns out that the optimization is also obtained for the same a. This 
can be seen by looking at a + xy 7 ^ and see that x = gives the maximum. 

Finally, if only one between §1 and § 2 is in (e, c), the other being smaller 
than e, we have to replace u_ )6 by t>_ j0 or u +>e by t> +j£ . In this case we then 
obtain the bound exp ( — |(e//?) 3 / 2 t + W(e 2 )) instead of the one in (15.901) . 
For < e <C 1 and t ^> 1, we then have the bound 

exp (-!(£/*:) 3/2 t). (5.91) 

Now, since e < §1 + s 2 < 2c, then £ > (si + s 2 ) 3 / 2 (/t/2c) 3 / 2 and -y/si + s 2 > 
t^^^/e/K. Thus we get 

(15311) < e -(si+ S2 ) 3 / 2 (e/2 C ) 3 / 2 /3 < e -(s 1 +s 2 ) e 2 i 1 /3 / ( 3c 3/2 K i/2 ) < e _( Sl+S2 ) , g 

for t large enough. □ 

A remark to point (2) of Proposition 15.41 The kernel is identically equal 
to zero as soon as the position we look is smaller than the initial position of 
the particle, thus the maximal value that the Sj's can take without the kernel 
being identically equal to zero is of order t. 
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Figure 7: Shape of Re(/ a;S ) along the real axis. The shown divergences are 
at —1/p, —1 and 0. 

Lemma 5.10. Let a > 0, < s < ft, and consider the function f a ^ s defined 
in ( 15. 85\) . F a ^ s has two stationary points at — 1 < < — 1/d < u+ }S < 0. 
Moreover, U- )S is strictly decreasing in s and u +)S is strictly increasing in s. 

Proof. In the limit case s = we have two solutions in (—1,0) but not strictly 
away from — 1/d, see (15.871) . First consider Let us take s > so > 0. By 
definition /4, So (w+, So ) = and lim «-o- fL,s ( u ) = ~oo, therefore f' aySo {u) < 
for u G {u +iSo , 0). On the other hand, 

f'a,s(u+,s ) = fa, so {u + ,s ) + ^ ~ > 0, S < S . (5.93) 

Therefore u +jS G (u +tSQ , 0), thus u +)S is strictly increasing in s. The result for 
m_ )S follows in a similar way. The shape of Re(/ a]S ) for a > 0, < s < ft is 
shown in Figure [71 

□ 

Lemma 5.11. Let a > and Si < ft. Then the following asymptotic expo- 
nential behaviors hold 

h(a, si) ~ exp (t/ a , 5l (tt-, s )), /i(a,Si) ^ exp ( - tf a>S2 (u +)S )). (5.94) 

Proo/. We will show that T_i = {u = -1 + («_,, + l)e'*,0 G [— 7T,7r)} and 
r = {v = —u +iS e ld ,6 G [— 7r,7r)} are steep descent paths^ of the integral 
defining I\ and I2 with maximal value at 6 = 0. The Lemma follows then 
from the fact that the paths r_i and r , which pass through u_ s and u+ s , 
are of length of order 1 in t. 

1 For an integral I = J dze l ^ z \ we say that 7 is a steep descent path if (1) Re(/(z)) 
is maximal at some z £ 7: Re(/(z)) < Re(/(zo)) for z e 7 \ {z } and (2) Re(/(z)) is 
monotone along 7 \ {zq} except, if 7 is closed, at a single point where Re(/) is minimal. 
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Let us start with J^a, s). Denote g{6) = Re(/ ajS (u)) for u G r_ 1; and 
r = m_ s + 1. Explicitly, 

g(B) = I In ((1 - p) 2 + p 2 r 2 + 2p(l -p)rcos(0)) 
Ol 

+- In (1 + r 2 - 2r cos(0)) + Constn (5.95) 
for some Constn independent of 9. Then 

A ' M - r MO) ( , ° m - . !,„ , m ) . (5.96) 



d0 v 7 Vl + r 2 -2rcos(0) ±_£ + ^r 2 + 2r cos(0) 

9 = corresponds to u — u- jS , which is a simple zero of d ^"^ , thus a saddle 
point, u = U- jS is a local minimum of Re(/(tt)) along the real direction, 
thus it is a local maximum along the path r_i at 9 = and for |0| <C 1, 
Re(/(u)) ~ — c0 2 for some c > 0. This implies that the term inside the 
brackets is strictly negative at 9 — 0. 

To see that T_i is a steep descent path, we have to check that the term 
inside the brackets is negative, as already shown for = 0. Thus we have to 
show that 

— ^ + — ^r 2 + 2r cos(0) < a' 1 (l + r 2 - 2r cos(0)) . (5.97) 
p 1 — p 

This rewrites as 

^(l+a-xl-lzf + r'g-j^). (5.98) 

This inequality holds for all because it holds for 9 = where the cos(0) is 
maximal. 

Now, consider I 2 (a,s). Let g{9) = Re(— f a , s ( v )) f° r ^ £ To an d denote 
r = — w+ )S . Then 

1 3 — s 

g{9) = -- ln(l + p 2 r 2 - 2pr cos(0)) + — — In(l + r 2 - 2r cos(0)) + Constn 

(5.99) 

for some Consti2 independent of 0. Then 

d *C) - rsin(S) ( ' _ _i_ . (5.100) 



d0 \ 1 + r 2 — 2r cos(0) p 1 + pr 2 — 2r cos(0) / 

If the term in the parenthesis is negative for all 0, then r is a steep descent 
path with maximum at u+ tS . This condition writes 

3 — s 

1 - + r 2 (l - p{3 - s)) > 2r cos(0)(l - 3 + s). (5.101) 

P 
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Since 1-/3+S > l-(3 = l-(d-l)p/(d-p)+(d-l)a > l-(d-l)p/(d-p) > 
for p G (0, 1), it follows that the inequality holds for all 9 if it holds at 9 = 0. 
But there the inequality holds since w +jS is a local quadratic minimum along 
the real direction of g(u). □ 



Proof of the bound of the binomial term 

Proof of Proposition 15.51 We divide the analysis into two cases. 

1) \si — s 2 \ < t 1 / 6 . In this case the results follows from the estimate (I5.2ip 
with O(e) < 0(t _1 / 6 ). This, together with the inequality 

leads to the desired bound. 

2) \si — s 2 \ > t 1 ^. Denote a as in f)5.20p and Xi — x 2 — 1 = aa. a ~ t 2 ^ 3 3> 1 
for large t. Then \si — s 2 \ > t 1 ^ 6 corresponds to \a — d\ > Cit -1 / 6 with 
ci = did — l)/(2n(r 2 — r\)). The a-dependent term in the l.h.s. of (15.71) is 

a \ / aa 



d-l 



-aa 

ag(a) 



e a9[a \ (5.103) 



with 

g(a) = (a ln(a) - (a - 1) ln(a - I) - a \n(d/{d - 1))) . (5.104) 

In (15.1031) the correction to e a9<y °^ behaves like 1/y/a, thus a prefactor 
1/V r 2 — r i wm always be in Const i. We have 

^ = ln(a/(a-l))-ln(d/(d-l)), ~~t~~^ = 7 ~ 77- (5-105) 
aa acr a (a — 1) 

Thus = only at a = d and is increasing in a. Thus for a G [1, d], 

9(a) < g(d) + A(a), f\(a) = -— L_( Q - d) 2 . (5.106) 

Moreover, for a G [l,Qf c ], 

^(a)<^(d)+/2(a), / 2 («) = - 4 ^ 1 _ 1) («-d) 2 (5.107) 

where a c is defined by /^(aj = lim^oo (/(a): a c = d+2d(d— 1) ln(<i/(<i— 1)). 
For a > a c , we set /3(a) = (a — d)f 2 (a c )/( y a c — d). Then g(a) < g(d) + / 3 (a). 
With these bounds on (7(a) we then get: 

a) 1 < a < d: e ag ^ < e a9 ^ e _(S2_S1 ) 2 7( 4 0*2-ri))^ 

b) d < a < a c : e ag{a) < e a 9( d ) e - (s2 - Sl)2/(8(r2 - ri)) , 

c) a > a c : e as(a) < e °-g(<i) e -\s2-s 1 \t 1 / 3 Kd 2 (a c )/(2d(d-i)in(d/(d-i))) _ 

Finally, by (15.1021) the desired bound holds for t large enough. □ 
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A Kernel and trace-class 

Predholm's series point of view 

One way of looking at (12. 6p is by simply considering the Fredholm series 
det(l - X-#FiX»)^({ui,... I u m }xB.) 
= ^2 ~ / dx i"-/ dx n det(i^(M ift ,x fc ;M i; ,x«))i< fc 



Kn- 



Then it is not difficult to see that the function is absolutely sum- 
able/integrable, so that the series is well defined. The required bound is 
easily obtained if we conjugate the kernel with 

p(u k ,x) = (l + x 2 ) 2k (A.l) 

as follows. We use 

p{u in xi) 



det {K(u ik ,x k \Ui v xi))i< Kl < n = det I K{u ik ,x k \u iv x h . , 

V P{ U i k i X k) J \<k,l<n 

■ (A.2) 

and the Hadamard bound, that is, the absolute value of a determinant of a 
n x n matrix with entries of absolute value not exceeding 1 is bounded by 
n n l 2 . The details are like in the proof of Theorem 12.51 



Operator point of view 

The second point of view is to consider (12.61) as a Fredholm determinant 
of an operator. The Fredholm determinant is well defined for trace-class 
operators. What we then mean by (12.61) is that there exists a conjugate 
operator (i.e., leading to the same determinantal measure) which is trace- 
class on 7i = L 2 ({ui, . . . , u m } x R). 

Consider the conjugate operator K\ given by the kernel 

Ki{u-i, si; u 2 , s 2 ) = K^m, si; u 2 , s 2 ) P , U2,S2 . . (A.3) 

p(«i,si) 

The choice of p{uk,x) is not at all unique. We use the following one: 

p(u k ,x) = (l + x 2 ) 2k . (A.4) 

Proposition A.l. The operator K\ as defined above is trace-class on 
H = L 2 ({u u . . .,u m } x R). 
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Proof. We can write K\ = XXz=i PkK\Pi where P& is the projection from TC 
into its subspace TCk = {/ € TC\f(ui, x) = 0,i ^ k}. Then 

m 

||#i||i.*< (A.5) 

k,i=i 

Now, HPfci^iPiHi^ = where K x ' is the operator on L 2 (R) with kernel 

Ki(uk, •; M;, •). Thus to prove that K\ is trace-class on TC, it is enough to prove 
that, for k, I = 1, . . . ,m, K^ 1 is trace-class as operator on L 2 (R). This is 
proven in Lemma [A. 21 and Lemma [A. 3 1 □ 

Lemma A. 2. The operator with kernel 

L k Mv) = (e^-^ )A )(x, W )^^l [:l .> Sfc] l [?/ > Si] l K>U;] (A.6) 
is trace-class on L 2 (R). 

Proof. For Uk < Ui nothing has to be proven since L^i = 0. Consider Uk > Ui 
(i.e., k > I) and define the two operators 

A(x, g ) = (e5^ A )(x >g ) ( ( ^P (A.7) 

and 

^y) = (e^ )A )(^y) ( \^ ) J w . (A.8) 

Then 

L M = P Sk ABP Sl (A.9) 

with P a (x) = l[ a > ]. The P's are projectors, thus ||P a . ||oo = 1- From this it 
follows 

\\L k ,t\\i < \\Ah\\B\\ 2 . (A.10) 

It is simple to prove that \\A\\2 < oo and ||P||2 < oo, since the 2-norm is 
easily bounded using its integral kernel. 

I 2 < f drdz 1 p -2{x-zf/{u k - Ul ) i 1 + ^ 2 ) 2(fc+0 1 

l2_ i R 2 27T(Ufc-Ml) (1 + X 2)2(fe+Z) (1 + x 2)2(fe-0 

(A.ll) 

By changing the variable z = y + x, it is easy to see that 

A z \ -2(x-z) 2 /(u k -u t ) (! + z ) ( + } < n nnstiQ c a in) 
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Therefore 

, dX (1 + z 2)2(ife-0 ^ 2 COI1Stl3 (A,13) 

since k — I > 1. In the same way we can check that ||-B|| 2 < oo, thus proving 
the result of the Lemma. □ 

Lemma A. 3. The operator with kernel 

M k j(x,y) = (e- U;A V= A )(a; !2/ )^^l [x > Sfc] l [?/ > S;] (A.14) 
is trace-class on L 2 (R). 

Proof. First of all, let P a be the projector on [a, oo), and set 

Then, for — L < min{-u fc , it holds that 

M M = P Uk Q k ,iP Ul = P Uk P-LQ k ,iP-LP Ul (A.16) 

and, using ||P a |U = 1, we get 

llM^Hx < \\P- L Qk,iP-L\\i. (A.17) 

Instead of using the projectors, we can think of Q k ,i as an operator on 
7i L = L 2 ([—L, oo), dx). The Airy function can be expressed as 

Ai(x) = — f d^e-^e^-^e^h 3 -^ 2 ^ (A. 18) 

for any a > 0. Let us set 
<*(0 = K 3 -^ 2 + K-^) 2 C, P = \^ -{uk-mfa+Huk-mf (A.19) 

and define the measures d/j,(x) = (1 + x 2 )~ 2l dx, du(x) = (1 + x 2 ) 2k dx. With 
these notations we have, for / e "Wl, 



|(/, Qk,if)n L \ 

and 



d//(x) / dis(y)f(x)Q k! i(x,y)f(y) 

-L J-L 



Qk,i(x,v) = f[ d^-^e ta ^e^ + y^- r \ (A.21) 



(A.20) 
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where we set r = — {u k — u\ — a). By choosing a > u k — U\ we have r > 
and we can exchange the integrals. Let us set 

/oo poo 
-L J~L 

(A.22) 

Then 



\{f,Qk,if)n L \ < f- ! dte-^\ gi (Z,a)\ \g 2 (i,a)\ 

< ^ f die-^{\ gi {^a)\ 2 +\g 2 {^a)\ 2 ) (A.23) 

where we used a 2 + b 2 > lab. The first term can be evaluated as follows. 

/ d^\ gi (^a)\ 2 (A.24) 



(POO \ / POO 

J ^ d^x)f(x)e-^e-™\ (J d^j/j/ge^e-^ 



d/i(x) / My)f(y)f(x)e- T ^ / aee-« 

— L J-L Jr 

/»oo POO 

" f d^x) / d /U ( 2/ )7M/(^)e- T( ^ ) e-^- y)2/4CT . 
The second term is computed in essentially the same way. We obtain 

< \(LQ k ,if)n\ < -r^={f,Gf) n (A.25) 

4\/ IT (7 



with G the operator with kernel 

G(x, y) = e-^e-^) 2 /^) ( (1 + + y2)2fc + (1 + ^-2^ + y2) -2i) 

(A.26) 

Qfc,z is Hilbert-Schmidt on Hl (use (IA.21I) ). Thus Qk,i is a compact 
operator. For any orthonormal basis {4> n } n >i we have 

E e /3 
|(0n,Qfc,i0n)wJ < V(0„,G0 n ) Wi . (A.27) 

n>l v n>l 



By ()A.25j) and (1A.26I) G is a positive and continuous operator, thus we have 



/oo POO 
dxG(x,x)= / da;e- 2T "((l+x 2 ) 4/ +(l+a; 2 )- 4fc ) < oo. 

(A.28) 

This holds for any orthonormal basis, thus the Lemma is proved. □ 
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B Krawtchouk orthogonal polynomials 

Let K n (x,p,t) = 2-^i(— n, — x; — t; 1/p) be the Krawtchouk polynomials, see 
e.g. [9]. To obtain an integral representation we use their generating function 

^(^K n (z,p,T)C= (l-^c)"(l + C) T - z - (B.l) 



n=0 

From this it follows that 



t\... _. i r .A 1 -^) ( l + T - 



nJ I U :,P,T) = -f r ^ ^ (B.2) 

where T is an anticlockwise oriented simple loop with as the only pole 
inside r . 

In the case of particles starting from i/j = —d(j — 1) we obtain, for 
k = 0,...,N- 1, 

*f w = (!-,)•(, (j^)"^^ - v-i)k,p,t + v. 

(B.3) 

We use the two identities 

(x + l)K n (x,p,N) (B.4) 
= (N + l)p(K n (x + l,p,N + 1) — K n+1 (x + l,p,N+ 1)) 

and 

(N+l-x)K n (x,p,N) (B.5) 
— (N + 1)(1 — p) (K n (x,p, iV + 1) + -P-K n+1 (x,p, N + 1)) 

recursively and eventually obtain, setting T — t + d(N — 1), 

d(N-l) 

*%(z)=u T {z) J2 S^K^p.T) (B.6) 

1=0 

with 

ur(z) = (1-p) T ( T JL -Y ( T ) (B.7) 



1 — pj \z 

the standard weight for the Krawtchouk polynomials Kj(z,p,T), 
j — 0, . . . , T, and the (not square) matrix 5 with entries 

--te);AE( (d -; )i )(T-/--; ) )te)" i <■»> 
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To obtain a family of polynomials (z),k = 0, . . . ,N — 1}, with $^ 
of degree k, and which satisfy ^2 z>0 $k ( z )^f( z ) = $k,j, we need to be able 
to invert the matrix S restricted to the entries < i,j < N — 1. The 
orthogonality relation of the Krawtchouk polynomials is, for < n, m < T, 

^u T (z)K m (z,p,T)K n (z,p,T) = (jz-^j ) <W- (B.9) 



z=0 

Therefore 



KM = f>^,p,T)(^) {jz^jsTi- (B.10) 

In the inverse of the matrix S = [5ij]o<i,.j,<JV-:b the terms in front of 
the sum are easily accounted for. In fact, if a matrix S has entries Sij = 
f(i)Sijg(j), then its inverse has entries S^J = g~ l (i)S7j f~ l (i)- In our case, 

f(i) = P~ l and g(j) = , and S h] = E,>o ( ( ^) (T^)^" with 

v 

Case d = 2 

In this case, set T = t + 2{N - 1) and oo T {z) = (1 - p) T (i^;) 2 Q- Then 



and 



1 - p/ V -p / VV V i - * - A 

2(AT-1) 

^(z)=uj t (z) SuKfanT), (B.12) 

= f>^,p,T)(^) (jz^fsrj. (B.13) 

The computation of the inverse of S is quite involved, but at the end of the 
week one finds, with < i,j < N — 1, 

l -Zl\ ' J ^ flz£ V,_, y- (2l/-2-i-j- A\ (J + A\ jW 



(B.14) 
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with 






2 (j + A)(2iV-2-z-j-A) 



the last being the expression for < A < j — i — 1. 

To obtain the functions $^ we plug the solution of the inverse of the ma- 
trix in the expression above and perform the sums. First one exchanges the 
path integral and the finite sums. Then the sums can be simplified by extend- 
ing them to — oo or oo depending on the case. This can be done because the 
extra terms turn out to be analytic functions leading to a zero contribution 
after integration. In the end we obtain $^ = 1 and, for 1 < k < iV — 1, 
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